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Abstract. Modifications of existing clouds by the exhaust 
of ships are well-known but inadequately quantified impacts, 
which could contribute to climate change. The perturbation 
of a cloud layer by ship-generated aerosol changes the cloud 
reflectivity and is identified by long curves in satellite im- 
ages, known as ship tracks. As ship tracks indicate a pollu- 
tion of a very clean marine environment and also affect the 
radiation budget below and above the cloud, it is important 
to investigate their radiative and climatic effects. Satellite- 
data from MODIS on Terra are used to examine a scene from 
10 February 2003 where ship tracks were detected close to 
the North American West-Coast. The cloud optical and mi- 
crophysical properties are derived using a semi-analytical re- 
trieval technique combined with a look-up-table approach. 
An algorithm is presented to distinguish ship-track-pixels 
from the unperturbed cloud pixels in the scene and from this 
the optical properties of the former are compared to those 
of the latter. Within the ship tracks a significant change in 
the droplet number concentration, the effective radius and 
the optical thickness are found compared to the unaffected 
cloud. The resulting cloud properties are used to calculate 
the radiation budget below and above the cloud. Assuming 
a mean solar zenith angle of 63° for the selected scene, the 
mean solar surface radiation below the ship track is decreased 
by 43.2 Wm~? and the mean reflectance at top of atmosphere 
(TOA) is increased by 40.8 Wm~?. For the entire analyzed 
scene the ship emission decreases the solar radiation at the 
surface by 2.1 Wm~? and increases the backscattered solar 
radiation at TOA by 2.0 Wm~?, whereas no significant effect 
on thermal radiation was detected. 
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1 Introduction 


Emissions by ships significantly contribute to the total bud- 
get of anthropogenic emissions. The exhaust gas emissions 
from ships include CO2, NOx, SOx, CO, hydrocarbons, and 
particulate matter (Eyring et al., 2005a). Recently, satel- 
lite measurements of NO2 showed enhanced NO2 column 
amounts along major international shipping routes (Richter 
et al., 2004; Beirle et al., 2004). Ship emissions are re- 
leased into the marine boundary layer and change the chem- 
ical composition of the atmosphere (Lawrence et al., 1999; 
Kasibhatla et al., 2000; Davis et al., 2001; Endresen et al., 
2003) and climate (Capaldo et al., 1999; Endresen et al., 
2003). The average sulphur fuel content of today’s world- 
merchant shipping fleet is 2.4% that results in a large amount 
of SOz2 and particulate matter emissions (EPA, 2000). In- 
ternational maritime conventions do not regulate the sulphur 
content of fuel for shipping. However, there are regional and 
local emission legislations, for NOx but also for SO, and par- 
ticulate matter (e.g., Entec, 2002; World Bunkering, 2004). 
During the last 50 years, the world’s ocean-going fleet and 
the ships’ total fuel consumption have increased consider- 
ably. Further increase of fuel consumption and emissions 
is expected in the future related to an increase in economic 
growth and sea borne trade (Eyring et al., 2005a, b). Ship 
emissions have been recognized as a growing problem by 
both policy makers and scientists (Corbett, 2003). 

In addition to the impact on tropospheric chemistry, parti- 
cle emissions from ships also change the physical properties 
of low clouds. This is the so-called indirect aerosol-effect, 
which has been observed in satellite data in many studies 
(e.g., Conover, 1966; Twomey et al., 1968; Radke et al., 
1989). Natural aerosol in a maritime environment consists of 
a small amount of particles mainly composed of sea salt with 
diameters greater than 1 wm (O’Dowd et al., 1997, 2004) and 
sulphate particles with diameters less than 1 yum (McInnes et 
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al., 1996). The low number of particles results in small op- 
tical thickness (Christopher and Zhang, 2002) and a small 
amount of aerosol available for cloud condensation. This 
limited number of cloud condensation nuclei is reflected in 
larger droplets and a smaller droplet number concentration 
in low-level stratiform clouds over the ocean compared to 
continental clouds (Miles et al., 2000). If additional aerosols 
are injected, the changes of the aerosol concentration and 
amount result in a change in the droplet number concentra- 
tion within the cloud (Facchini et al., 1999), depending on 
the solubility and size of the injected aerosol particles. Par- 
ticles and their precursors from ship emissions are able to 
act as cloud condensation nuclei (CCN) in the water-vapour 
saturated environment of the maritime cloud or can change 
the surface tension due to the solubility. Which particulate 
matter of ship emissions is able to act as CCN is not entirely 
clarified, but especially the high sulphur content of the fuel 
may be an important factor for the modification of clouds, 
because the resulting SOx is able to act as CCN (Seinfeld 
and Pandis, 1997). Amount and size of these particles de- 
pends on the fuel and also the kind of combustion, but can 
possibly result in a higher droplet concentration (Twomey et 
al., 1968; Twomey, 1974) and consequently in a change of 
reflectivity of the maritime cloud. 

Measurements in the Monterey Ship Track Experiment 
confirm this hypothesis (Durkee et al., 2000a). In the plumes, 
Hobbs et al. (2000) observed an increased amount of those 
particles which can act as CCNs in a supersaturated envi- 
ronment as well as the related increase in the cloud droplet 
concentration and the decrease in the cloud effective radius. 
Slight changes in the particle size distribution were reported 
by Öström et al. (2000). On the other hand, there was no 
evidence for the influence of sea salt particles (Durkee et al., 
2000b) from the wake of the ship on the clouds and also no 
significant influence of heat and moisture from the ship stack 
(Hobbs et al., 2000). So, these measurements confirm the 
role of emitted particles as CCNs and especially the impor- 
tance of SOx. 

Modifications of clouds by aerosols result in a signifi- 
cantly higher droplet concentration. This leads to an in- 
creased scattering, resulting in larger reflectances. This in- 
creased reflectance however, might be partially reduced due 
to the presence of highly absorbing particles (e.g. soot) in 
the plume exhaust. The increased reflectivity is even higher 
in the near infrared, because here, the ratio of absorption to 
scattering is strongly depending on the droplet size (Coakley 
et al., 1987; Kokhanovsky et al., 2004a). 

In this study the modification of clouds and the influ- 
ence of the ship exhaust on the radiation budget of a given 
scene are examined. Satellite data are used to retrieve cloud 
properties and their modifications due to ship emissions. 
Section 2 gives an overview of the selected satellite scene. 
Section 3 explains the approach for extracting the polluted 
(ship track) pixels and the algorithms to analyze microphys- 
ical and optical properties of the cloud. They are based 
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on the semi-analytical cloud retrieval algorithm SACURA 
(Kokhanovsky et al., 2003) combined with look-up-tables for 
thin clouds calculated with the libRadtran radiative transfer 
package (Mayer and Kylling, 2005). The work concentrates 
on low clouds because only those are directly affected by 
ship emissions. The section also addresses how foam and 
turbulent waves caused in the ship wake may affect the re- 
flectances and hence the extraction of ship tracks. The last 
part of this section describes calculations to determine the 
impact on the radiation budget. 

The results of this study are presented in Sect. 4. The cloud 
mask separating ship-track-pixels and no-ship-track-pixels is 
shown in Sect. 4.1. Changes in the effective radius, opti- 
cal thickness and droplet number concentration due to the 
exhaust plume are presented in Sect. 4.2. The resulting op- 
tical parameters obtained from the image analysis are used 
in Sect. 4.3 to calculate the radiation below and above the 
cloud and to estimate the influence of the modified clouds on 
the radiation close to the surface and also at the top of the 
atmosphere (TOA). 


2 Selected ship track scene 


For our analysis a particular and adequate satellite scene 
from Terra-MODIS (King et al., 1995) was selected us- 
ing a resolution of 1km. The scene from 10th February 
2003, close to the West Coast of North America (153° W to 
120° W and 40°N to 60°N), exhibits a number of anoma- 
lous cloud lines in the stratiform clouds over the ocean. Fig- 
ure 1 shows the reflectance in MODIS-data for channel 2 
(0.85 um) and Fig. 2 the reflectance for channel 7 (2.13 um). 
These two channels are also used for the ship track detec- 
tion. As explained above, the curves are more pronounced 
in the 2.13 ~m-channel. The time of the satellite overpass is 
20:25 UTC (11:25 local time), the solar zenith angle ranges 
between 55° and 79°, depending on the location in the im- 
age. The observed type of clouds is typical for this region in 
the Northern Pacific (Norris, 1998a, b; Whang et al., 2000). 
The scene contains the Northern Pacific and on the east- 
ern edge, the American west coast is visible. Large ports 
like Vancouver and Seattle (see Fig. 1) are located on this 
coastline. A comparison with | x1 degree vessel traffic den- 
sity for the year 2000 (Endresen et al., 2003) for the area 
(Fig. 3) shows high density of ship traffic for the domain 
of the satellite scene. Any impact of air traffic is elimi- 
nated by taking into account the brightness temperature of 
the clouds, using the MODIS atmospheric window chan- 
nel 31 (11 um, Fig. 4). The temperatures indicate low warm 
clouds for the anomalous cloud lines and therefore the impact 
of aircraft emissions is negligible. Meteorological conditions 
taken from ECMWF reanalysis data (ERA-40-project, Sim- 
mons and Gibson, 2001) for a geopotential height of 700 m 
above the ocean show a high pressure field over the northern 
Pacific with an anticyclonic wind field (Fig. 5). The curves 
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Fig. 2. Same as Fig. 1, but for 2.13 um-channel. 


are consistent with ship-tracks being formed by ship-traffic 
coming from or going to the great harbours of the West coast 
of North America being turned by the wind field. In conclu- 
sion it is very unlikely that this anomalous cloud lines have a 
source other than ship traffic emissions. 
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3 Methods 


An automated algorithm to detect ship tracks in this scene has 
been developed in order to study how ship emissions modify 
clouds. Figure 6 shows a flow diagram of the ship track algo- 
rithm. The algorithm distinguishes between ship-track-pixels 
and no-ship-track-pixels, that is, polluted and unpolluted pix- 
els. In addition, optical and physical properties of the low 
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Fig. 3. Vessel traffic density of the total fleet for the year 2000 
between 140° E to 120° W and 15° N to 75° N. Data are taken from 
Endresen et al. (2003). 


clouds are calculated and the derived optical properties are 
used to estimate the influence on the radiation budget. 


3.1 Ship track mask algorithm 


The upper part of the flow diagram in Fig. 6 shows the cri- 
teria that are used to identify ship-track-pixels. Sea pixels 
are first separated from land pixels by means of a land-sea- 
mask. Then a filter is applied to extract only those cloudy 
pixels, which are susceptible for the exhaust of ships. To de- 
termine a threshold for cloudy pixels, a histogram analysis 
of the reflectance was performed. For a pixel to be clas- 
sified as cloudy, its reflectance at 0.65 um (Channel 1) has 
to be higher than Rthreshhoia=0.1. For a pixel to be classi- 
fied as low cloud, a temperature threshold was used. Dur- 
kee et al. (2000b) present measurements of the height of low 
stratiform clouds which are modified by ship emissions. The 
cloud top height never exceeded 1500m in these measure- 
ments. Taking into account the temperature of the ocean and 
a lapse rate for a wet-adiabatic environment like the maritime 
boundary layer, the cloud top height can be estimated from 
the brightness temperature in channel 31 (11 um). To sepa- 
rate low cloud pixels from high cloud pixels the temperature 
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difference AT between the sea surface temperature (T ocean) 
and the temperature of the cloudy pixels had to be less than 
10.5 K. The sea surface temperature Tocean was derived from 
nearby cloud free pixels, while AT has been calculated as- 
suming a wet-adiabatic lapse rate of 0.7 K/100 m for the mar- 
itime boundary layer and a maximum top altitude of 1500 m. 
The temperature criterion also excludes the detection of air- 
craft contrails, as the algorithm only detects warm clouds. 
The remaining cloud covered pixels were classified as low- 
cloud-pixels. 

Finally, pattern recognition, in particular Canny-Edge- 
filtering, was used in order to identify changes in the cloud 
reflectance, which indicate ship tracks. Enhancement of the 
edges was achieved by dividing the values of the satellite data 
from channel 2 (0.85 um) and channel 7 (2.13 um). Chan- 
nel 7 is a good indicator for cloud modifications and therefore 
ship tracks, because the reflectivity strongly depends on the 
ratio of absorption to scattering in the near infrared, which 
can be seen in Fig. 2. Channel 2 helps to increase contrast 
for pixels with increased scattering. Using both channels, 
the contrast for pixels with a significant influence on the ra- 
diation, is increased. 

The first step, Canny-Edge-filtering, does not distinguish 
between ship-track-pixels and natural cloud edges. To fur- 
ther select the ship track relevant pixels automatically the 
following criteria were applied in addition: 


— Width of ship track: The number of pixels between two 
edge-pixels should not exceed a threshold, chosen for 
narrow curves; 


— Ratio of inside/outside-pixel: The length and width of 
connected pixel-clusters should not exceed a threshold; 


— Ratio of length/width: The number of neighbouring and 
connected pixels in an area should not exceed a thresh- 
old. 


The thresholds were optimized manually to minimize the 
false detection of clouds, by excluding edges that were obvi- 
ously not ship tracks but natural clouds, e.g. for this scene we 
found an optimum at a width of 10 pixels between the edges. 
The result of the ship track mask algorithm described above 
is first a differentiation of the low-cloud-pixels from high 
cloud-pixels and second a separation of low-cloud-pixels into 
the two categories ship-track-pixels and no-ship-track-pixels. 

In Sect. 4 the sensitivity of the different parameters of the 
cloud mask are described. To minimize the false detection of 
ship track-pixels, the algorithm was optimized empirically 
to identify ship-track-pixels in the scenes, which contribute 
significantly to the radiative forcing. These are defined as 
radiatively relevant ship-track-pixels. 


3.2 Impact of uncertainties due to surface reflectance 


Ships create wakes marked by foam with increased re- 
flectivity. For thin, semi-transparent clouds this increased 
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Fig. 4. Temperature of the selected scene in Kelvin, derived from satellite data (11 um). 


reflectance might be misinterpreted as a modification of the 
cloud parameters. To analyse reflectance changes produced 
by the wake of the ship, sensitivity simulations have been 
carried out with libRadtran (Mayer and Kylling, 2005) at 
0.85 um, 1.64 um, and 2.13 wm. 


The sensitivity simulations were calculated using a solar 
zenith angle of 60°, a satellite zenith angle of 10°, a relative 
azimuth angle of 45° and cloud optical properties represen- 
tative for the conditions of the examined satellite image (see 
Sect. 4.2.). A first calculation (SIM1) was designed to simu- 
late a “normal” cloud over an almost black Lambertian ocean 
surface (spectral reflectances of 0.03 at 0.85 um and 0.0 at 
1.64 um and 2.13 um), using an effective radius of 13.2 um 
and a optical thickness of 20.7 for a no-ship-track pixel. A 
second calculation (SIM2) was designed to simulate an hypo- 
thetical ship wake filling the whole pixel, and used an ocean 
surface with spectral reflectances of 0.6 for 0.85 um and 0.3 
for 1.64 um and 0.1 for 2.13 um (Kokhanovsky, 2004b). In 
SIM2 the same cloud optical parameters as in SIM1 have 
been used. The third calculation (SIM3) is similar to SIM1, 
but uses cloud optical parameters typical for a ship track 
rather than for an unpolluted cloud: In SIM3 an effective 
radius of 10.1 um and a cloud optical thickness of 34.6 were 
used. 


Figure 7a shows the calculated differences in reflection 
between SIM1 and SIM2 and betweeen SIMI and SIM3. 
The change in reflectivity in the visible wavelength channel 
(0.85 um) is almost the same for changes in surface albedo. 
In this case, changes in reflectivity due to ship tracks could 
not be distinguished from changes in surface reflectance, e.g. 
due to foam in a wake. However, for the longer wavelength 
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Fig. 5. Wind field from ECMWF reanalysis data on 10 February 
2003 at 18:00 UTC. 


channels (1.64 um and 2.13 um) the change in reflectivity 
between ship track and foam is obvious. 

Adding the information of a wavelength in the near in- 
frared reduces the possibility to interpret changes in surface 
albedo due to foam as changes in cloud optical parameters. 
Calculations at 2.13 um showed the best signal-to-noise ratio 
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Fig. 6. Flow chart of the ship-track-algorithm that has been used in 
this study to identify ship tracks and calculate microphysical, opti- 
cal and radiative properties. 


such that the use of this channel is particulary helpful for ship 
track detection. 

The same three sensitivity calculations were repeated for a 
thinner cloud with a cloud optical thickness of 5.0 and an ef- 
fective radius of 13.2 um for “normal” optical cloud parame- 
ters and a cloud optical thickness of 10.0 and an effective ra- 
dius of 10.1 for a ship track cloud. This is shown with same 
albedo assumptions as before in Fig. 7b. Here, the ocean 
albedo has a remaining influence also on 2.13 um that could 
be misinterpreted as a change in cloud contrast by the ship- 
track-mask-algorithm. 


To exclude possible misdetection, in addition a smaller 
scene has been selected (hereafter: “smaller selected scene”), 
where the cloud optical thickness of the unpolluted clouds 
near the ship tracks is very high. The optical thickness in 
this subset is 10 to 30 in the no-ship-track-clouds. Therefore, 
effects of changes in contrast due to foam or turbulence for 
2.13 um are in any case negligible for the selected area. Also 
the amount of cloud edges and fractional cloud cover is re- 
duced in this smaller selected scene. Therefore, the detection 
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of ship tracks in the smaller selected scene is more reliable 
than the entire scene. 


3.3 Cloud properties retrieval 


We developed our own algorithm in order to be able to ap- 
ply the retrieval to different satellite instruments. This new 
algorithm combines the semi-analytical cloud retrieval algo- 
rithm SACURA (Kokhanovsky et al., 2003) and look-up- 
tables (LUTs) for thin clouds calculated with the libRadtran 
radiative transfer package (Mayer and Kylling, 2005). An 
advantage of the algorithm is that it can be applied to dif- 
ferent satellite instruments with channels in the near infrared 
(e.g. Terra-MODIS, AATSR). 

The optical and microphysical parameters of the cloud 
were derived from MODIS channel 2 (0.85 um) and chan- 
nel 6 (1.64 um). 1.64 um was selected because the smaller 
absorption of liquid water enables more accurate results for 
the SACURA-retrieval. The two cloud retrieval algorithms — 
SACURA and LUTs — are different, but both derive the cloud 
optical thickness and also the effective radius reff, defined by 
the ratio of the third to second moment of the particle size 
distribution and therefore indicating a change of the ratio of 
volume to surface in the particle size distribution. 

SACURA is a semi-analytical approach based on the 
fact that the probability of photon absorption by droplets 
in the visible and near-infrared spectral regions is low 
(Kokhanovsky et al., 2003). The technique relies upon 
asymptotic equations of the radiative transfer theory valid 
for the values of the single scattering albedo wọ close to 
one, which allows a much faster calculation of reflectances 
than a numerical solution of the radiative transfer equation. 
SACURA was compared to other cloud retrieval algorithms 
(Nauss et al., 2005) and the retrieved values agreed well. 
SACURA calculates the effective radius from the observed 
reflectances and uses this value to derive optical thickness 
and liquid water path, assuming a vertically homogeneous 
cloud. The ocean surface reflectance for SACURA was as- 
sumed to be 0.03 for 0.85 um and 0.0 for 1.64 um. The ad- 
vantages of SACURA are its speed and high accuracy and 
also an interpolation is not necessary. A disadvantage of 
SACURA is the increasing error for decreasing cloud opti- 
cal thickness. 

Therefore, for clouds with an optical thickness of 5 and be- 
low, a look-up-table approach and an interpolation is used, to 
have a wider range of applicability in terms of the cloud op- 
tical thickness. The look-up-table (LUT) is calculated using 
the radiative transfer code libRadtran (Mayer and Kylling, 
2005), assuming mid-latitude standard atmosphere and a 
cloud top height of 1000m. The scattering properties of the 
cloud are based on Mie-calculations and the assumed surface 
reflectance was the same as for SACURA. Reflectances were 
calculated for a wide range of optical thicknesses and effec- 
tive radii. A best-fit method is used to find the cloud optical 
thickness and rege for the observed reflectances. 
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Fig. 7. Changes in simulated reflectivity (TOA) due to changes in surface reflection by foam or changes of cloud optical properties. Panel (a) 
shows a thick cloud (no ship track: 20.7 for optical thickness and 13.2 wm for effective radius, ship track : 34.6 and 10.1 um), panel (b) 
shows a thin cloud (no ship track: 5.0 and 13.2 um, ship track: 10.0 and 10.1 um). 


The columnar droplet number concentration for both re- 
trievals is calculated via effective radius and cloud optical 
thickness by assuming a gamma droplet size distribution with 
a coefficient of variance of 0.37 as proposed by Deirmend- 
jian (1969). Calculations of droplets per volume (N) were 
performed using a hypothetical vertical homogeneous cloud 
of a thickness 500 m, which is a reasonable value for low ma- 
rine stratiform clouds (Cahalan et al., 1994; Han et al., 1998). 
This calculation does not influence the relative changes and 
resulting values per volume are useful for comparing with 
measured cloud droplet concentrations, e.g. Miles (2000). 

Threshold criterion for calculations with either SACURA 
or LUT was the reflectance of channel 2 (0.85 um), using 
a value, which approximately fits to the cloud optical thick- 
ness around 5. The cloud optical thickness is calculated at 
0.55 um for the LUT and 0.85 um for SACURA, but cloud 
optical thickness is approximately wavelength-independent 
in this region. Calculations with both retrievals around 
+/—0.1 of the reflectance-selection criteria near the thresh- 
old showed a correlation factor of 0.73 for the optical thick- 
ness with a mean relative deviation of 11.3%, and 0.91 for 
the effective radius, with a mean relative deviation of 6.7%. 
The differences are within the uncertainty of the theoretical 
model, e.g. the assumption of an unpolluted vertically and 
horizontally homogeneous cloud within a pixel. 


3.4 Estimating the impact on the radiation field 


The derived optical parameters were used to estimate 
changes in solar radiation for the areas below and above the 
cloud as well as the thermal outgoing radiation by radia- 
tive transfer calculations. Optical thickness and the effec- 
tive radius have been applied to create look-up-tables for the 
solar flux via the radiative transfer code libRadtran (Mayer 
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and Kylling, 2005), by using the built in k-distribution by 
Kato et al. (1999) to calculate integrated solar irradiance with 
the solver disort2 (Stamnes et al., 1988) for the wavelength 
range of 0.24 um to 4.6 um. The down-welling irradiance 
at the surface and the up-welling flux at TOA were calcu- 
lated for the mid-latitude winter atmosphere (AFGL-TR-86- 
0110, Anderson et al., 1986). The different distributions of 
cloud optical properties were considered by using these look- 
up-tables to calculate the solar flux for every pixel and tak- 
ing into account the local solar zenith angle. The cloud top 
height was chosen to be 1000 m and the cloud-bottom height 
was 500m. The optical properties of the clouds were cal- 
culated according to Mie theory. The mean values for all 
low-cloud-pixels, ship-track-pixels and no-track-pixels were 
determined (results see Sect. 4.3), to estimate the impact of 
ship tracks on both, the solar radiation at the surface and the 
backscattered radiation at TOA. 

In a similar manner to the calculations described above, 
changes in thermal radiation resulting from changes in cloud 
optical properties were calculated (0.2 um to 100 um), using 
the correlated-k distribution by Fu and Liou (1992). Estima- 
tion of exact cloud top height is uncertain for these clouds, 
but possible influence on cloud top height due to aerosol was 
included. The thermal radiation of every pixel was calcu- 
lated by taking into account the actual cloud top height. This 
quantity was determined by matching the measured bright- 
ness temperatures at 11 um (Fig. 11) to an appropriate tem- 
perature profile by assuming a lapse rate of 0.7 K/100 m. 
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Fig. 8. Results of the ship track mask, plotted over the satellite data shown in Figs. 1 and 2. Ship-track-pixels are indicated in blue. 


Table 1. Sensitivity analysis: Detected ship-track-pixels for variations of the cloud-mask-parameters, compared to the 6.7% detected ship- 


track-pixels that are used for the examination. 


Change of parameter (%) 


Detected ship-track-pixels (%) 


Canny-Edge Parameter —10/10 5.3/8.6 
—50/50 2.2/31.5 
Width of ship track —10/10 6.2/7.1 
—50/50 4.0/8.2 
Ratio inside/outside-pixel —10/10 6.6/6.8 
—50/50 6.6./6.8 
Ratio of length/width —10/10 6.6/6.8 
—50/50 6.1/7.1 


4 Results 
4.1 Detected ship tracks 


Figure 8 shows the resulting analysis of the particular scene 
derived from the ship track mask algorithm described in 
Sect. 3.1. The resulting ship-track-pixels are shown in blue 
together with the satellite data of Fig. 1 and 2. A small num- 
ber of pixels marked as ship-track-pixels are not inside a ship 
track and almost a third of the ship tracks seem not to be de- 
tected. However, the ship tracks visible at the shorter wave- 
lengths are almost all detected and those are the relevant pix- 
els for the solar radiative budget. By using the combination 
of channels 7 and 2 only those ship-track-pixels are detected, 
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which significantly change the backscattered radiation at top 
of atmosphere. For the whole scene, 6.7% of the low-cloud- 
pixels were identified as ship-track-pixels and 93.3% as no- 
ship-track-pixels. For the smaller selected scene 8.6% ship- 
track-pixels and 91.4% no-ship-track-pixels were identified 
by the algorithm. 


To test the sensitivity of the algorithm on the cloud mask 
algorithm parameters, we have varied the parameters for 
Canny-Edge filtering, the width of the ship track, the ratio 
of inside and outside pixels as well as the ratio of length to 
width (see Sect. 3.1) by +/— 10% and 50%. Table 1 shows 
the detected ship-track-pixels in percent of the entire scene 
for the different input parameters. Results are most sensi- 
tive to changes in the Canny-Edge parameter. A change of 


www.atmos-chem-phys.net/6/4925/2006/ 


M. Schreier et al.: Impact of ship emissions on properties of marine stratus 


4933 


Table 2. Mean values of cloud parameters and standard deviation (STD) in parenthesis. 


low-cloud-pixels 


no-ship-track-pixels — ship-track-pixels 


All low-cloud-pixel 


effective radius (um) 13.0 (4.0) 
optical thickness 21.4 (11.1) 
droplet number (cm7?) 85 (61.7) 

Smaller selected scene 

effective radius (um) 11.6 (3.9) 
optical thickness 24.4 (10.5) 
droplet number (cm73) 115 (60) 


13.2 (4.0) 10.1 (2.5) 
20.7 (10.4) 34.6 (14.5) 
79 (82.8) 210 (112.6) 
11.9 (3.8) 8.9 (2.3) 
23.8 (10.2) 30.9 (11.8) 
101 (46) 255 (101) 


Table 3. Changes in mean values of retrieved cloud optical properties by changes in cloud-mask-parameters. 


ship track pixels noship track pixels 

Change eff. opt. N eff. opt. N 
of radius thick. (cm7) radius thick. (cm7?) 
param (%) (um) (um) 

Canny-Edge Parameter 10 10.3 34.2 198 13.0 20.9 81 
50 11.2 28.9 140 12.8 20.8 81 
Width of ship track 10 10.0 34.4 215 13.1 20.8 80 
50 9.5 33.5 225 13.4 21.0 79 
Ratio of inside/outside 10 10.0 34.5 212 13.1 20.7 80 
50 10.0 34.8 214 13.1 20.7 80 
Ratio of length/width 10 10.0 34.5 212 13.1 20.7 80 
50 10.0 34.6 213 13.1 20.7 80 


+/—10% results in a loss of 1.4%/increase of 1.9%, respec- 
tively. Results are also sensitive to the width of the ship track 
defined as the pixels from one edge to another. A value of 10 
pixels was used as default in the algorithm. Changing this to 
9/5 pixels results in a decrease of 0.5%/4%, respectively. In- 
creasing the value to 11 or 15 pixels results in an increase of 
0.4% or 1.5%. Results are not sensitive on the ratio of inside 
and outside pixels and the ratio of length to width. 


4.2 Microphysical and optical cloud properties 


The cloud retrieval algorithm (Sect. 3.3) was used to cal- 
culate optical and microphysical parameters of low clouds 
(Fig. 9). Ship emissions result in the decrease of the effec- 
tive radius (Fig. 9a). The increase of the optical thickness is 
also significant (Fig. 9b). Obvious is a high increase in the 
droplet concentration (Fig. 9c) compared with the unpolluted 
cloud scene. 

Figure 10 shows the smaller selected area to allow a more 
detailed analysis of the ship track characteristics. In this area 
no cloud-edges are included and also the solar zenith angle 
(SZA) varies between 60° and 66° so can be assumed as al- 
most constant. In addition, the optical thickness of the no- 
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ship-track clouds is mostly higher than 10 (see Sect. 3.2) and 
therefore changes in reflectivity due to foam can be excluded. 
A significant decrease of the effective radius from 12 um to 
6 um is visible across the ship-track-pixels (Fig. 10a). The 
optical thickness of unpolluted clouds is about 20 to 30 and 
is increasing in the track up to 45 and higher (Fig. 10b). Also 
the change in the droplet number concentration from around 
100cm~? up to 800 cm7? is substantial (Fig. 10c). 


Table 2 summarizes the mean values of the various pa- 
rameters for the low-cloud-pixels, both for the entire area 
as well as for the selected area for ship-track-pixels and no- 
ship-track-pixels. The decrease in the effective radius from 
13.2 to 10.1 um for the entire area and from 11.9 to 8.9 um 
for the smaller area is evident and also an increase in cloud 
optical thickness from 20.7 up to 34.6 for the entire area and 
from 23.8 up to 30.9 for the selected area is observed. 


Table 3 provides an analysis of retrieved mean optical 
cloud properties (effective radius, optical depth and number 
concentration), if the cloud mask parameters are varied ac- 
cording to Table 1. The mean optical parameters are only 
sensitive to the edge filtering (Canny-Edge-parameter) and 
the width of the ship track. 


Atmos. Chem. Phys., 6, 4925-4942, 2006 


4934 


M. Schreier et al.: Impact of ship emissions on properties of marine stratus 


150.0 
100.0 
50.0 


g- 


aS 


a) 30.0 b) 80.0 
28.0 75.0 
26.0 70.0 
1500 24.0 1500 65.0 
22.0 60.0 
20.0 veel 
Cc 
18.0 ON E 
a 45.0 8 
= 1000 16.0% — 1000 too È 
x aS ws O: 
ee 14.05 `å = 
ia 2° 35.0 $ 
42, = 
Ss 30.0 & 
10.0 25.0 9 
500 8.0 500 20.0 
6.0 15.0 
4.0 10.0 
20 5.0 
o E 0.0 o 0.0 
o 200 400 600 800 1000 o 200 400 600 800 1000 
pixel pixel 
c) 800.0 
750.0 
700.0 
1500 650.0 Q 
fe} 
600.0 V. 
a 
550.0 5 
500.0 3 
o 
450.0 ? 
_— 1000 Q 
g 400.0 9 
a. O 
350.0 8 
300.0 6 
250.0 S 
500 200.0 2 


(0) 


0 200 400 


600 
pixel 


0.0 


800 1000 


Fig. 9. Effective radius (a), cloud optical thickness (b) and concentration of particles (c) (cm~3) derived from the MODIS channels 2 and 6 


for the analysed scene. 


To evaluate our results presented above, we compare the 
parameters cloud optical thickness and effective radius with 
the same parameters in freely available MODIS cloud stan- 
dard product (King et al., 2003; Platnick et al., 2003). For the 
entire scene an effective radius of 13.2 um for the no-ship- 
track-pixels and 10.6 wm for the ship-track-pixels is calcu- 
lated for the MODIS standard product. The optical thick- 
ness increases from 17.8 in no-ship-track pixels to 29.8 in 
ship-tracks-pixels. In the smaller selected scene, the effec- 
tive radius decreases from no-ship-track-pixels (11.7 um) to 
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ship-track-pixels (8.6 um) and the optical thickness increases 
from 17.6 to 23.3. Overall, the MODIS cloud standard prod- 
uct gives very similar results for the effective radius, whereas 
the optical thickness in our retrieval is slightly larger, which 
is consistent with Nauss et al. (2005). 


Figure 11 shows the resulting statistical distributions of the 
parameters comparing no-ship-track-pixels and ship-track- 
pixels for all low-cloud-pixels (A, B, C) and for the smaller 
selected area (a, b, c). If all low-cloud-pixels are consid- 
ered (Fig. 11A), a change of the maximum of the effective 
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Fig. 10. Effective radius (a), cloud optical thickness (b), and droplet number concentration (c) derived from the MODIS channels 2 and 6 


for the smaller selected area. 


Table 4. Mean values of the transmitted solar flux at the surface and reflected radiation at the top of atmosphere, STD in parenthesis. 


low-cloud-pixels 


no-ship-track-pixels 


ship-track-pixels difference 


Elow—cloud Eno-ship—track Eship—track Ejow—cloud—Eno—ship—track 
All low-cloud-pixel 
Surface (Wm) 127.6 (61.1) 130.1 (61.0) 78.7 (36.1) —2.5 
TOA (Wm?) 336.2 (55.9) 335.4 (56.9) 353.5 (26.2) 0.9 
Smaller selected scene 
Surface (Wm~?) 113.4 (28.3) 115.5 (27.5) 92.6 (27.6) —2.1 
TOA (Wm?) 349.8 (16.2) 348.2 (15.6) 365.1 (13.3) 1.6 


radius towards smaller radii can be seen, whereas the distri- 
bution is changed towards higher values for the cloud optical 
thickness (Fig. 11B). The maximum of the droplet number 
distribution curve shows a significant shift of the maximum 
from 70 cm7? to 160 cm~? (Fig. 11C). Droplet numbers are 
higher than 500cm~°, in contrast to the ones outside of the 
tracks. The distribution of the cloud optical thickness inside 
the tracks shows values again above 50 (Fig. 11b), which do 
not appear in the distribution of the outside pixels. The ship- 
track-pixels of the selected area also show a narrower distri- 
bution for the effective radius and the cloud optical thickness 
confirmed by a decrease of the standard deviation. 
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In summary, ship emissions of aerosols and their precur- 
sors result in an increase of the droplet number concentration 
and a decrease in the effective radius in the cloud. A larger 
number of smaller droplets for a fixed liquid water content 
results in the observed increase in the cloud optical thick- 
ness. The results from satellite data analysis agree with pre- 
vious studies, e.g. from the MAST-experiment (Ostrom et al., 
2000; Hobbs et al., 2000; Ferek et al., 2000). 
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Fig. 11. Statistical distributions of the parameters for the no-track-pixel and the ship-track pixel for the entire scene (capital letters) and the 


smaller selected scene (lowercase letters): 


(A, a) effective radius, (B, b) cloud optical thickness, (C, ¢) droplet number concentration. 


4.3 Radiative changes at surface and TOA 


4.3.1 Radiative effects 

Figure 12 shows down-welling solar radiation at the surface 
below the clouds and reflected solar radiation at TOA, calcu- 
lated by the method described in Sect. 3.4. In this calculation 
the optical parameters derived as described in Sect. 3.3. and 
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depicted in Fig. 9, as well as the local solar zenith angle is 
taken into account for each pixel. As the local solar zenith 
angle varies over the scene from approx. 50° to 80°, the sim- 
ulated surface and TOA irradiance are influenced by both, 
the SZA and the optical parameters. Below a ship track the 
solar radiation is significantly reduced compared to outside 
the track because the reflectance of the cloud in the track is 
higher. The reflected solar radiation at TOA above the ship 
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Fig. 12. Down-welling solar radiation at the surface below the clouds (a, in Wm?) and (b) reflected solar radiation on TOA due to the 


clouds (b, in Wm7~2). 


Table 5a. Mean values of solar radiation at the surface and TOA for a fixed solar zenith angle of 63°, standard deviation in parenthesis. 


fixed sza=63° low-cloud-pixels 


Ejow—cloud 


no-ship-track-pixels 


Eno— ship—track 


ship-track-pixels difference 


Eship—track Elow—cloud-Eno—ship—track 


All low-cloud-pixel 


Surface (Wm?) 124.2 (45.0) 126.3 (44.8) 
TOA (Wm72) 332.1 (38.2) 330.1 (37.8) 
Smaller selected scene 

Surface (Wm72) 113.5 (44.1) 115.5 (44.2) 
TOA (Wm) 343.5 (37.4) 341.3 (37.3) 


83.1 (28.9) —2.1 
370.9 (23.5) 2.0 
93.6 (27.9) —2.0 
364.9 (22.8) 2.2 


tracks is significantly enhanced. Table 4 shows the mean val- 
ues of the down-welling solar radiation at the surface and the 
reflected solar radiation at TOA, averaged over all ship-track- 
pixels (Eship—track), no-ship-track-pixels (Eno—ship—track), and 
low-cloud-pixels (Ejow—clouds) for the whole scene and the 
smaller selected area. The resulting transmitted energy be- 
low the ship-track-pixels is reduced by 2.1 and 2.5Wm7?, 
for the smaller selected area and the entire scene respectively. 
The corresponding changes at TOA are an enhancement of 
the reflected solar radiation by 1.6 Wm~? for the smaller se- 
lected area and 0.9 Wm? for all low-cloud-pixels. 


When absolute radiation values (ship-track-pixels vs. no- 
ship-track-pixels) are compared, caution is required since 
the incident solar radiation at TOA depends on geographic 
position (i.e. on solar zenith angle) and ship tracks are not 
uniformly distributed over the selected scene. To quantify 
the magnitude of the solar zenith angle impact on the so- 
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lar flux, the variation of the average incoming irradiance at 
TOA was simulated with the radiative transfer model for all 
pixels, all ship-track- and all no-ship-track-pixels for both, 
the full scene and the smaller selected area. The difference 
between all low cloud pixels and the no-ship-track-pixels 
is about —0.6Wm~? (all low-cloud-pixel) to —1.6 Wm~? 
(smaller selected area). In addition to that, the cloud re- 
flectance, transmittance as well as atmospheric absorption 
also depend on solar zenith angles. In summary, the impact 
of the solar zenith angle variation over the scene is of simi- 
lar order of magnitude as the combined effect of changes in 
optical properties and solar zenith angle effects. 


To separate the impact of changes in cloud parameters 
from ship tracks on the radiation field without the uncertain- 
ties represented by the solar zenith angle variations, radiative 
transfer calculations were performed assuming a mean solar 
zenith angle of 63° (Table 5a) for all pixels. The calculated 
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Table 5b. Same as Table 5a, but for a solar zenith angle of 50°. 


fixed sza=50° low-cloud-pixels 


Elow—cloud 


no-ship-track-pixels 
Eno- ship—track 


ship-track-pixels difference 


Eghip—track Elow—cloud-Eno— ship—track 


All low-cloud-pixel 


Surface (Wm7~) 213.9 (77.7) 217.5 (77.3) 
TOA (Wm72) 442.4 (65.2) 439.0 (64.6) 
Smaller selected scene 

Surface (Wm?) 195.3 (76.2) 198.9 (76.3) 
TOA (Wm72) 461.4 (64.0) 457.7 (63.8) 


142.7 (49.8) —3.6 
508.1 (40.0) 3.4 
160.9 (48.1) —3.5 
497.5 (39.1) 3.7 


Table 5c. Same as Table 5a, but for a solar zenith angle of 70°. 


fixed sza=70° low-cloud-pixels 


Elow—cloud 


no-ship-track-pixels 
Eno— ship—track 


ship-track-pixels difference 


Eship—track Elow—cloud-Eno— ship—track 


All low-cloud-pixel 


Surface (Wm?) 80.6 (29.1) 82.0 (28.9) 
TOA (Wm?) 258.6 (24.8) 257.3 (24.6) 
Smaller selected scene 

Surface (Wm7~) 73.8 (28.5) 75.1 (28.6) 
TOA (Wm?) 266.1 (24.4) 264.7 (24.3) 


54.0 (18.7) a 
284.0 (18.1) 13 
60.9 (18.1) 3 
280.3 (14.8) 1.4 


values in Table 5 now only depend on the cloud optical prop- 
erties, as the influence of a varying solar zenith angle has 
been eliminated. 


Assuming a constant solar zenith angle, the mean value 
of the down-welling solar radiation at the surface below the 
ship-track-pixels is reduced by 43.2 Wm~? due to the change 
in the cloud parameters, compared to the no-ship-track-pixels 
in the entire scene. This change is comparable to the calcu- 
lated discrepancy of 63 Wm~* between maritime and conti- 
nental clouds as proposed by Kiehl (1994) for a similar solar 
zenith angle. At TOA, ship-track-pixels reflect 40.8 Wm~ 
more than the no-ship-track-pixels for the entire scene. For 
the selected area the reduction of the down-welling solar ra- 
diation at the surface due to ship tracks is 22.1 Wm~? and the 
increase in TOA is 23.6 Wm~°. 


The net radiative effect of the change in cloud properties 
due to ships for the particular scene is estimated by calcu- 
lating the difference of the absolute radiation values between 
all low-cloud-pixels (Ejow—cloud) and the no-ship-track-pixels 
(Eno—ship—track). According to these values, the solar radia- 
tion at the surface is reduced on average by 2.1 Wm~? by the 
ship emissions and additional 2.0 Wm? are reflected back 
at TOA. Table 5a also shows the values for the smaller se- 
lected area with an on average decrease in the radiation of 
2.0 Wm“? at the surface and an increase in the outgoing ra- 
diation at TOA of 2.2 Wm~? due to the ship emissions. 
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For the same spatial distribution of cloud and ship track 
properties, the solar radiation has been also calculated for 
other solar zenith angles to get an estimate of the depen- 
dence of the radiation changes by ship tracks on the solar 
zenith angle (Tables 5b—d). For a smaller solar zenith angle 
of 50° the scene average radiation decreases to —3.6 Wm~? 
at the surface and increases to 3.4 Wm? at TOA. For larger 
mean solar zenith angle of 70° the scene average change is 
—1.4Wm~? at the surface and 1.3 Wm~? for the backscat- 
tered radiation at TOA. For an even large solar zenith angle 
of 80° the change is only —0.5 Wm~? for the surface and 
only 0.5 Wm~? for TOA. In summary, smaller SZA result in 
an amplification of the solar radiation changes on the surface 
and TOA, whereas for high solar zenith angles the effect will 
vanish. Ship tracks can have al long lifetime, e.g. Durkee et 
al. (2000b) give an average estimate of 7.5 h but they also re- 
port a high variability. Therefore, the solar zenith angle vari- 
ation has to be included when estimating the radiative impact 
of ship tracks. A possible change in thermal radiation was 
also analysed taking into account the changes in cloud opti- 
cal properties and also changes in cloud top height. A small 
reduction of 0.7 Wm~? at TOA was found by comparison of 
non-ship-track-pixels and ship-track-pixels. This results in a 
negligible change of less then 0.1 Wm~? for the entire scene. 
The results are similar for the smaller selected scene. 
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Table 5d. Same as Table 5a, but for a solar zenith angle of 80°. 
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fixed sza=80° low-cloud-pixels 


Elow—cloud 


no-ship-track-pixels 
Eno- ship—track 


ship-track-pixels difference 


Eghip—track Elow—cloud-Eno— ship—track 


All low-cloud-pixel 


Surface (Wm~7) 30.4 (10.9) 30.9 (10.8) 
TOA (Wm72) 134.6 (9.4) 134.1 (9.3) 
Smaller selected scene 

Surface (Wm7~) 27.9 (10.7) 28.4 (10.7) 
TOA (Wm72) 137.5 (9.2) 136.9 (9.2) 


20.4 (7.1) —0.5 
144.3 (5.8) 0.5 
23.0 (6.8) —0.5 
142.9 (5.6) 0.6 


The results show that anthropogenic emissions from ships 
influence maritime clouds and affect the radiative energy 
budget at the surface and also the reflected radiation at TOA. 
Ship tracks show an increased reflection of solar radiation 
compared to unaffected clouds. For the scene investigated 
in this work, the net-effect is a reduction by 1 to 3 Wm as 
depending on how the solar zenith angle impact is estimated. 

Not included in this analysis is the diurnal cycle of mar- 
itime stratiform clouds (Cahalan et al., 1994; Norris, 1998b) 
and the lifetime of ship tracks, also influencing the radiative 
behaviour. 


4.4 Impact of cloud mask parameters on radiation values 


In addition to the sensitivity of the cloud mask parameters 
on the number of detected ship-track-pixels (Sect. 4.1.) and 
microphysical properties (Sect. 4.2.), here we estimate the 
sensitivity of the cloud mask parameters on changes in radi- 
ation. 

As can be seen in Fig. 11, the cloud optical parameters 
are not equally distributed which results in a non-linear cor- 
relation of the resulting loss or augmentation of pixel as de- 
scribed in Table 1. Table 6 shows the sensitivity of the re- 
sulting loss of radiation at TOA. Changing the parameters of 
the Canny-edge-filtering by 10% results in relatively small 
changes in the radiation. But for a variation by 50%, results 
change significantly. The edge filtering picks up only the 
most reflective pixel and the reduced pixel amount of 2.2% 
is only the highly radiative ship-track-pixel, corresponding 
to a value of 7.5 Wm. On the other side, the increased de- 
tection of 31.5% can be explained by a higher detection of 
previously no-ship-track-pixels as ship-track-pixels resulting 
ina decrease of 0.7 Wm~*. Therefore, a high amount of ship- 
track-pixels is required for such calculations. For only 2.2% 
of the ship-track-pixel of the entire scene, the calculation of 
the net effect of the scene would increase the estimation of 
the net effect by 5.5 Wm. For the remaining factors of the 
algorithm, variation works vice versa: a reduction of the de- 
tected pixel results in a reduction of the radiation value. In 
this case, the edge filtering did already pick up the radiatively 
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relevant pixels. The remaining factors can only reduce the 
mean value by reducing the pixel amount. For these factors, 
the value has only an important change, if the width of the 
ship track is reduced or increased by 50%. 

In conclusion, the cloud mask algorithm extracts the most 
radiatively relevant ship-track pixel via the edge filtering. 


5 Summary and conclusion 


An algorithm has been developed to determine ship tracks 
from satellite data in an automated way. A scene on 10 
February 2003 was chosen to extract the optical and mi- 
crophysical cloud modifications from ship emissions using 
Terra-MODIS satellite data. A combination of the semi- 
analytical approach SACURA with a look-up-table for op- 
tically thin clouds was used to calculate cloud properties. 

It was demonstrated that the wake of a ship might influ- 
ence the ship track detection if changes of the surface would 
reach pixel size. Adding a channel in the near infrared helps 
to exclude errors by foam. 

It has also been shown that ship emissions modify existing 
clouds by decreasing the effective radius, while they increase 
droplet concentration and optical thickness. On average, the 
optical thickness was increased from 20.7 up to 34.6 and 
the effective radius was decreased from 13.2 um to 10.1 um. 
The calculated average droplet number concentration in- 
creased from 79 up to 210 cm7?. The results agree with mean 
values calculated by MODIS standard cloud product and also 
with the theory and experiment (Öström et al., 2000; Hobbs 
et al., 2000): Low clouds of the maritime boundary layer 
have less cloud condensation nuclei than clouds over land; 
in consequence, this results in larger droplet radii for simi- 
lar water content and dispersion of droplet size distributions. 
Injection of aerosols and their precursors by ships results in 
more CCNs causing the mean droplet radius to decrease and 
the droplet number concentration to increase. 

The derived parameters were used to calculate changes in 
the radiative energy budget below and above the cloud. The 
mean values show a decrease of 43.2 Wm~? for the surface 
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Table 6. Change of results for entire scene, if variation of the cloud mask parameters is done as in Table 1. 


Change of parameter (%) 


backscattered Radiation (Wm?) 


Canny-Edge Parameter —10/10 
—50/50 
Width of ship track —10/10 
—50/50 
Ratio inside/outside-pixel —10/10 
—50/50 
Ratio of length/width —10/10 
—50/50 


2.4/1.6 
7.5/0.7 
1.8/2.1 
1.1/2.4 
1.9/2.0 
1.7/2.1 
1.9/2.0 
1.8/2.1 


radiation below ship tracks and an increase of 40.8 Wm~? for 
the increased reflectivity at TOA. If the whole low-cloud area 
with 6.7% ship-track-pixels is taken into account, a decrease 
in the radiation at the surface of 2.1 Wm~? and an increase of 
2.0 Wm~ in backscattered solar radiation was found, when 
assuming a constant solar zenith angle of 63 ° for the scene. 
A small reduction of thermal radiation was found at TOA 
with almost no effect for the entire scene. 

Changes in radiation have also been calculated for a 
smaller area of the selected scene. The results show sim- 
ilar changes, supporting the radiation effects calculated 
for the entire scene: a decreased radiation at the sur- 
face of 2.0 Wm”, an increased backscattered radiation of 
2.2Wm?. 

The results showed that modifications of clouds by inter- 
national shipping can be an important contributor to climate 
on a local scale. Further studies are needed to assess the 
global impact of ship-track formation on climate. 
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